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Abstract 

A realistic computer model for polymorphic crystallization (i.e., initial and final phases with identical compositions), which 
includes time-dependent nucleation and cluster-size-dependent growth rates, is developed and tested by fits to experimental data. 
Model calculations are used to assess the validity of two of the more common approaches for the analysis of crystallization data. 
The effects of particle size on transformation kinetics, important for the crystallization of many systems of limited dimension 
including thin films, fine powders, and nanoparticles, are examined. © 1997 Elsevier Science S.A. 

Keywords: Computer modeling; Crystallization kinetics; Microstructure 


1. Introduction 

An ability to predict and to control the crystalliza- 
tion of metallic and silicate glasses is often critical for 
the preparation and maintenance of desirable mi- 
crostructures. Further, the quantitative analysis of ki- 
netic crystallization data can provide information about 
the magnitude and temperature dependence of nucle- 
ation and growth rates, essential parameters for materi- 
als design. Kinetic data are generally obtained by 
monitoring some physical parameter that changes dur- 
ing crystallization, such as the electrical conductivity. 
The interpretation of such data, however, is often un- 
clear; the electrical resistivity, for example, can be ex- 
tremely sensitive to the transformation microstructure 
[1]. The latent heat evolved during crystallization scales 
linearly with the rate of transformation, making calori- 
metric measurements such as differential scanning 
calorimetry (DSC) and differential thermal analysis 
(DTA) attractive. These techniques are sensitive, rela- 
tively easy to make and require small sample quantities. 
In many cases, however, experimental studies of crystal- 
lization suffer from a lack of methods for quantitative 
analysis. This is particularly true for nonisothermal 
data, where most methods are based on erroneous 
assumptions. Even for isothermal studies, however, 
proper account is frequently not taken of the transfor- 
mation microstructure, making many conclusions sus- 
pect. 

Analytical methods for the analysis of isothermal and 
nonisothermal studies of glass crystallization have been 


criticized on theoretical and experimental grounds [2- 
4]. Experimental [5,6] and theoretical [7] studies of the 
Kissinger method [8], the most commonly used tech- 
nique for the analysis of nonisothermal data, have 
demonstrated that it can be used with some confidence 
for transformations involving growth only. Within this 
short space it is not possible to review and critique the 
many analysis schemes that exist; a partial review' of 
these can be found in [2]. Here, a brief review is 
provided of our studies of polymorphic crystallization 
of silicate glasses, where the initial and final phases 
have the same composition. The most common analyti- 
cal methods of analysis, the Kissinger and Ozawa [9] 
methods, are examined using this model. It should be 
emphasized that a silicate glass is chosen only because 
nucleation and growth rates and required thermody- 
namic parameters are known; the conclusions drawn 
also apply to crystallization of metallic glasses. 

Accurate data for the time-dependent volume frac- 
tion crystallinity are assumed. Experimental details 
such as proper baseline corrections for calorimetric 
data are not discussed. Since time-dependent nucleation 
behavior is central to the model development, this is 
first reviewed briefly in Section 2. 

2. Non-steady-state nucleation rates 

Generally, the nucleation rate of crystals in a glass is 
taken to be constant at the steady-state value, 7 s . If 
true, the number of nuclei produced in a unit volume, 
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Time (s) Cluster Size, n 

/u'f 'k ^ Nun ’ bcr of nuc,el as a function of time, illustrating non-steady-state nuclcation rates and showing the definition of the induction time 
(b) schematic illustration of cluster evolution in the kinetic model for nuclcation. 


N vi should then scale linearly with time. In many 
silicate and some metallic glasses, however, a nonlinear 
behavior is often found (Fig. 1(a)), which indicates an 
initially slow nucleation rate that approaches 7 s after 
long annealing times. Such time-dependent, or tran- 
sient, homogeneous nucleation has been observed in a 
wide variety of condensed phases [10]; examples of 
heterogeneous transient nucleation have also been re- 
ported in the crystallization of silicate [11] and metallic 
[12] glasses. Studies of transient nucleation can provide 
insight into the dynamics of cluster evolution underly- 
ing nucleation phenomenon. Further, transient nucle- 
ation can determine phase formation and stability [7, 1 3] 
and often influences the final microstructure. 

Time-dependent nucleation can be understood within 
the classical theory of nucleation [14] where clusters of 
the new phase grow or shrink by the addition or loss of 
a single entity (monomer) at a time (Fig. 1(b)), follow- 
ing a series of bimolecular rate reactions. This cluster 
evolution is described by a system of coupled differen- 
tial equations of the form: 

<iN„. 

fa ~ i i — (N„,A ■„ + N„,k + ) + N n t l k, l . , 

( 1 ) 

where and k, ~ are the forward and backward rate 
constants, respectively, which are proportional to the 
diffusion coefficient in the initial phase, and N„, is the 
number of clusters at time t containing n monomers. 


The rates and cluster densities are governed by the 
work of cluster formation, W n , which for spherical 
clusters containing n monomers and having a sharp 
interface with the initial phase can be written approxi- 
mately as 

W n = nSfi 4- (36/r) 1 Y 2 V V (2) 

where S/.i is the Gibbs free energy per monomer of the 
new phase less that of the initial phase, v is the molecu- 
lar volume and a is the interfacial energy per unit area. 
The nucleation rate, /„,, is defined as the flux in cluster- 
size space; in the most general case, it is a function of 
both time and the cluster size at which it is measured, 

kij kt nt Ic n N n + 1 jk „ + 1 ( 3 ) 

Given sufficient time, the cluster density evolves by 
the mechanism described in Eq. (1) and approaches a 
constant value, the steady-state distribution, (V*,, pro- 
ducing a constant, steady-state rate. Since many metal- 
lic and silicate glasses are formed by rapid quenching 
from the melt, there is often insufficient time to main- 
tain this steady-state distribution of clusters, however, 
resulting in a final cluster distribution that is more 
similar to that at a higher temperature, having a low'er 
cluster density at the critical size and therefore a de- 
pressed nucleation rate. With time, the distribution 
evolves to the steady-state one, giving rise to the char- 
acteristic time-dependent nucleation rate. 
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3. Johnson -Mehl - A vr ami - Kolmogorov (JMAK) 
analysis 


In the very early stages of polymorphic crystalliza- 
tion (i.e., no compositional change) of glasses, overlap 
between the different transformed regions can be ig- 
nored and the volume fraction crystallized can be writ- 
ten as 




'4 f 

Jo LJr 


8(0 At' 


dr 


(4) 


Here, V u is the sample volume, /(r) is the time-depen- 
dent nucleation rate and g(t') is the time-dependent 
growth rate. As the crystallized regions become larger, 
they eventually impinge and alter the transformation 
kinetics. If the sample size is much greater than any 
individual transformed region and if growth proceeds 
homogeneously throughout the sample, with spatially 
random nucleation, then the effect of impingement on 
the kinetics can be computed statistically to yield the 
familiar JMAK expression [15,16], i.e., x(t) = 1 — 
exp( — .v c (r)). For interface-limited growth and a con- 
stant nucleation rate: 


x(t) = 1 — exp( — 4 /t/£ V/3 V u ) (5) 

Often, a general form of this type is assumed to de- 
scribe isothermal phase transformations in other cases: 

x(t)= 1 - expo (A' 0") 

The parameters describing the reaction kinetics, such as 
the nucleation and growth rates, are now contained 
within an effective kinetic parameter, A, while the expo- 
nent, iu called the Avrami exponent, provides some 
information about the dimensionality of the tiansfor- 
mation, i.e., whether it is one-, two- or three-dimen- 
sional and whether it is interface-limited or 
diffusion-limited. A simple manipulation of Eq. (6) 
shows that for isothermal cases, these values are deter- 
mined by plotting the measured data, x(t) 9 as 


In 


In 


1 


l 2C (/ ) 


= n ln(A) + n ln(/) 


(7) 


A plot of ln[ln(l - x) ~ l ] versus ln(0 should then give a 
straight line with slope ft and intercept // ln(A). 

Continuously curved lines in Avrami plots have led 
to attached importance for local values’ of the Avrami 
exponent [17]. Generally, however, such behavior sig- 
nals a violation of a fundamental assumption of the 
JMAK theory. A familiar example in crystals is grain 
boundary nucleation which violates the assumption of 
spatially random nucleation [18], for example. Exam- 
ples of similar violations in glass devitrification include 
surface crystallization [1] and an inhomogeneous distri- 
bution of quenched in nuclei [19]. It is important to 
note that when applied in situations consistent with the 
original assumptions, the JMAK analysis is valid. 


4. Finite sample size effects 


Samples for DSC/DTA studies are often small, in the 
form of powders, thin films or ribbons. For rapid 
nucleation in reasonably large samples, possible finite- 
size effects on the transformation kinetics can be ig- 
nored. For low nucleation rates, however, or extremely 
small sample sizes, as for non-consolidated nanoparti- 
cles, the size of an individual transforming region is of 
order the size of the sample, thus violating a fundamen- 
tal assumption of the JMAK analysis. Possible effects 
on the crystallization kinetics have been considered 
only recently [20-22]. 

The probability for nucleation in small samples is 
readily obtained using Poisson statistics; growth pre- 
sents a more difficult problem. Given a small ensemble 
of crystallizing regions and assuming isotropic growth, 
some crystallites will intersect the sample surface prior 
to complete transformation and growth will cease, thus 
violating the assumption of unimpeded growth that led 
to the JMAK equation. A recent study [21] shows that 
if the regions of significant nucleation and growth do 
not overlap (valid for many silicate glasses) the finite- 
size-corrected volume fraction transformed from N nu- 
clei can be written approximately in terms of a scaled 
time, k , 


x(t) = 1 - exp 


, 9 „ m(N) , 

n[ — +—H 2 *('> 


( 8 ) 


where K = gt/R, g being the growth velocity of the 
crystal region, R the sample particle radius, and Mhc 
time; m(N) is equal to the larger of the two quantities 
(5.26 0.26 N) or 1. Particle shape effects are easily 
included, using an effective value for k c = (VJV 0 ) 1 , 
where V c is the extended volume of a growing crystallite 
and V 0 is the volume of the sample. Predictions from 
Eq. (8) agree with computer simulations of the transfor- 
mation of a finite ensemble of noninteracting spins on a 
lattice (Fig. 2(a)- (c)), while those from the JMAK 
expression compare poorly. The Avrami exponent ob- 
tained by a fit of the JMAK equation to the finite-size- 
corrected data (Fig. 2(d)) is initially low, due to the 
inapplicability of the JMAK analysis, and appioaches 
the expected value of w = 3 as the number of clusteis 
becomes large. 


5. Realistic computer model for isothermal and 
nonisothermal crystallization 

A computer model of polymorphic crystallization 
that includes time-dependent nucleation, cluster-size-de- 
pendent growth rates and finite-size effects is developed 
and applied to the crystallization of lithium disilicate 
glass. 
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0.5 1.0 1.5 2.0 

Reduced Time (k) 



Fig. 2. The fraction transformed as a function of reduced time, k, for (a) one, (b) two and (c) six crystallites in a particle ((a) (c) taken from pi]) 
The simulation results were obtained on a spherical lattice of radius 50 units. Eq. (8) is used for the finite-size calculation, which provides the best 
lit to the lattice simulation, (d) The Avrami exponent computed from a JMAK analysis of the finite-size-corrected data. 


5. 1. Development of computer model 

The transformation time is divided into a series of 
short isothermal intervals of duration At , over which 
new nuclei appear and previously generated ones grow. 
The coupled differential equations (i.e., Eq. (1)) are 
solved using a finite difference method [23], in which 
each isothermal interval is divided into a large number 
of smaller intervals, St, and the number of clusters of 
size n at the end of each interval is computed by 


N„, t + st — N„ y/ + St 


d/ 


( 9 ) 


where dN„Jdt is given by Eq. (1). Given the 
time-dependent nucleation rate is computed directly 
using Eq. (3). By this method, no assumptions of 
thermal history are made. At the end of each At 
interval, nuclei generated in previous intervals are 
grown assuming a size-dependent growth rate [13]: 


g(R) = C 


1 6D 



13 


sinh 


2k T 



( 10 ) 


Here, k is Boltzmann’s constant, / is the jump distance, 
and AG V is the change in Gibbs free energy per volume 
of the new phase less that of the initial phase. For 
nonisothermal transformations, the duration of the 
isothermal interval. At, is inversely related to the scan 
rate, S. The extended volume fraction transformed, x c , 
is calculated at the end of each interval by summing 
over all previous m intervals: 


n-TT I 



( 11 ) 


Finite-size effects are included as discussed in Section 4. 
Since for many silicate glasses, the nucleation and 
growth are well separated, all crystallites are approxi- 
mately the same size, allowing the scaled time, x c , to be 
obtained directly from the extended volume fraction for 
an infinite system, \ c jl), i.e., a~ c = (v e _, f/)/A r j , ? . Surface 
crystallization is frequently important for silicate and 
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Fig. 3. The effect of scanning through the nucleation zone at different rates on 

and (□); calculation ( ). For 850-1170 particles: data (•); calculation (- - 

assume ellipsoidal particles with eccentricity (1, { ) (from [25]). 


the DSC/DTA peak parameters. For 425 450 particles: data (■) 
). Calculations in (a) and (c) are for spherical particles; (d) (e) 


metallic glasses. For the glass studied here, surface 
nucleation is sufficiently rapid that the nucleation step 
can be ignored [24]; the rate of transformation due to a 
spherical growth front originating on the surface can 
then be assumed. Taking the growth velocity for the 
shell as constant and equal to that predicted from Eq. 
(10) as r -» x , g(ao), the volume fraction transformed 
by surface growth alone is 


A’ 



fl- g(x)/ 
R 


( 12 ) 


Calculating the volume transformed by surface growth, 
and taking account of the overlap between surface, x s , 
and volume, x v , crystallization, the time-dependent vol- 
ume transformed is given by 


x(t) = x s + x v (l - -Y s ) V*) 

Assuming that the rate of change of the enthalpy, 
dHjdt scales linearly with the rate of volume trans- 
formed, the DSC/DTA signal can then be computed as 
a function of time and/or temperature as 


DSC Signal: 


x{T i + ST)-x{T i ) 


(14) 


5.2. Application to crystallization of lithium disilicate 
glass 

Measured parameters were assumed to model the 
crystallization of lithium disilicate glass. The fitting 
parameters were held constant at the values detei mined 
from direct nucleation and growth data; they were not 
optimized to produce better agreement with the experi- 
mental data for crystallization. DSC/DTA peak 
parameters, such as peak temperature, T v> and peak 
height, 7 P , are strong functions of particle size [24] and 
sample history; the computer model presented here is 
capable of taking both effects into account. As an 
example, Fig. 3 compares the model predictions with 
experimental data for peak parameters for two sizes of 
spherical glass particles that were first scanned through 
the region in which the nucleation rate is significant 
(400-500°C) at rates, 0> n , between 2 and 15°C nun 1 , 
and were subsequently scanned through the crystalliza- 
tion peak at a rate, ® c , of 15°C min" 1 . Symbols are 
used for the data; the lines are the calculations. Results 
from two independent DSC/DTA measurements of 
these glasses are presented, providing a measure of the 


St 
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Fig. 4. (a) and (b): Calculated relative peak heights, / p , for glasses that had been scanned through the nucleation zone ( ) or held isothermally 

for a corresponding amount of time at the peak temperature of the steady-state rate ( ) compared with experimental data ((•) scanned 

glasses; ( ) isothermally treated glasses). Calculations in (a) were made assuming transient nucleation; steady-state rates were assumed for (b). 
Corresponding calculated time-dependent nucleation rates are shown in (c) and total number of nuclei in (d). The solid and dashed lines are for 
samples scanned through the nucleation zone at 4> n = 4 and 0.5°C min — 1 , respectively; the dotted and dot-dot-dash lines are for samples annealed 
at the peak nucleation temperature for 18.75 and 150 min, respectively. The steady-state nucleation rate is also shown and the nucleation and 
crystallization zones are indicated. 


experimental error. Though the predicted trends for 
spherical particles (Fig. 3(a)-(c)) follow the data, 
slightly improved agreement is obtained if ellipsoidal 
particles are assumed (Fig. 3(d)— (f)). Optical mi- 
croscopy studies of the particles show that this assump- 
tion is reasonable. 

Time-dependent nucleation can significantly enhance 
glass formation [13] and glass stability [7], due to the 
smaller number of quenched-in nuclei. As is illustrated 
in Fig. 4, however, the inherited non-steady-state clus- 
ter distribution can also result in unexpected transfor- 
mation behavior. Here the computed and 
experimentally measured DSC/DTA peak heights are 
compared for glasses that were either first scanned 
through the nucleation zone (400 500°C), or annealed 
isothermally at the peak temperature for an amount of 
time equal to that required to scan over the range of 
significant nucleation (425-500°C). Both sets of glasses 
were subsequently scanned through the DSC/DTA 


peak at 15°C min -1 . The measured peak heights are 
approximately equal, but fall faster than predicted from 
steady-state calculations (Fig. 4(b)), indicating a 
stronger dependence on the number of nuclei produced 
than expected from the temperature dependence of the 
steady-state rate alone. The experimental and calcu- 
lated results are almost indistinguishable from the time- 
dependent nucleation calculation (Fig. 4(a)). 

Earlier calculations by us predicted that if transient 
nucleation effects were important, the nucleation rate 
could rise above the steady-state nucleation rate during 
nonisothermal devitrification, peaking at a temperature 
higher than the peak of the steady-state rate [7], A 
similar behavior is expected here following both 
isothermal and nonisothermal anneals (Fig. 4(c)); two 
peaks are even predicted when <F n = 0.5°C min - h Fig. 
4(d) shows the number of nuclei produced as a function 
of temperature for different annealing treatments; this 
complex behavior is a result of the relaxation of the 
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cluster distribution inherited from the quench. From 
steady-state arguments, the number of nuclei should 
scale linearly with the annealing time in the nucleation 
zone; the predicted time-dependent behavior is more 
complicated. For example, for O n = 0.5°C min 1 , most 
of the nucleation occurs in the nucleation zone. Though 
approximately the same number of nuclei are produced 
for the corresponding case of a 150 min anneal at the 
peak temperature of the steady-state rate, these are 
nucleated at a temperature above the upper limit of the 
nucleation zone. Anneals within the nucleation zone 
then directly produce nuclei, but equally important, 
they result in differing degrees of relaxation of the 
cluster distribution toward the steady state one. The 
observed agreement between the isothermal and non- 
isothermal annealing treatments is therefore only fortu- 
itous. It is, however, a dramatic demonstration of the 
often subtle effects of time-dependent nucleation on 
phase formation and stability. This also further demon- 
strates that simple interpretations of DSC/DTA data, 
even in apparently simple situations, are often not 
possible. 

5.3. Critique of Kissinger and Ozawa analysis 
techniques 


Having established the validity of the model, it can 
be used to make meaningful evaluations of commonly 
used methods of analysis for DSC/DTA data. Peak 
profiles for DSC/DTA data were first computed using 
known kinetic and thermodynamic parameters for 
lithium disilicate glass. These data were then analyzed 
in the usual way and the derived parameters were 
compared with those used for the calculations. Two 
method of analysis for DSC/DTA data are most com- 
mon. The Kissinger method [8], 



where Q is the DSC scan rate and T p is the peak 
temperature (K), is often used to estimate the effective 
activation energy 7 of the transformation, E. The Ozawa 
method [9] 


d(log[ln(l - a)]) 
d log Q 


(16) 


is used to estimate the Avrami coefficient, n, which can 
provide information about the transformation mecha- 
nism. The validity of these expressions was checked by 
application to numerically generated DSC/DTA trans- 
formation peaks for computed, quenched, glasses using 
known thermodynamic and kinetic parameters. DSC/ 
DTA data were generated for infinite samples and for 
cases where surface and finite-size effects were included. 


In all cases, straight lines were obtained for the 
Kissinger and Ozawa plots. As illustrated in Fig. 5(a), 
the activation energies obtained from Kissinger plots of 
the computed DSC/DTA data (given by the slopes of 
the lines) were a strong function of the sample thermal 
history and whether surface and sample size effects 
were included. Large values of n were obtained from an 
Ozawa analysis of the more rapidly quenched glasses, 
indicating that time-dependent nucleation is dominant 
during the transformation. A value of n near 3 for the 
glasses quenched at 0.0 1°C s -1 suggests that the trans- 
formation proceeds primarily by growth on quenched- 
in nuclei. This w 7 as verified by noting that the activation 
energies computed from a Kissinger analysis become 
more similar (Fig. 5(a)), and the Avrami coefficients 
approach 3 for all glasses that were first saturated with 
nuclei by annealing in the nucleation zone (400 
500°C). 

Since nucleation and growth are widely separated in 
temperature for lithium disilicate glass, the activation 
energies obtained from a Kissinger analysis should be 
related to those of the growth rate, particularly if it is 
described by an Arrhenius temperature dependence. 
For the more realistic non-Arrhenius growth rates in 
this glass, however, the situation is more complicated. 
As illustrated in Fig. 5(b), showing the local activation 
energy for growth as a function of temperature, activa- 
tion energies obtained by a Kissinger analysis corre- 
spond roughly to the average value of the activation 
energy over the range of transformation temperatures 
for the different scan rates, though the value is not 
located precisely at the midpoint of the range. As 
illustrated, changes in the activation energy resulting 
from the introduction of finite-size effects or from 
thermal annealing are due to a change in the tempera- 
ture range of transformation. Samples of finite-size are 
transformed at lower temperatures due to a more rapid 
transformation of the finite volume while low tempera- 
ture annealing simply results in more nuclei that collec- 
tively transform the sample earlier. Even for this 
relatively simple case of polymorphic crystallization 
with separation between nucleation and growth, a 
proper interpretation of activation energies derived 
from a Kissinger analysis is difficult to make, particu- 
larly when information is desired about the influence of 
nucleation from anneals below the DSC/DTA transfor- 
mation temperature. For cases where the nucleation 
and growth regions overlap considerably, a Kissinger 
analysis gives meaningless results. Since the nucleation 
rate has a strongly non-Arrhenius temperature depen- 
dence, a primary assumption of the derivation is vio- 
lated. Because the Ozawa method avoids this 
assumption, the estimated Avrami coefficient more ac- 
curately reflects the importance of nucleation at tem- 
peratures below the transition temperature in the most 
rapidly quenched glasses. 
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Fig. 5. (a) Kissinger analysis for as-quenched glasses and glasses that held been annealed in the nucleation zone (400 500°C); (b) local activation 
energy for growth velocity. 


Most techniques for the analysis of nonisothermal 
crystallization data rely on similar assumptions to those 
made for the Kissinger analysis and therefore have the 
same level of validity. How then can this data be better 
analyzed? As discussed in Section 5.2, it should be 
possible to distinguish between crystallization mecha- 
nisms based on a qualitative analysis of the peak 
shapes. An iterative computer fit to the peak parame- 
ters, using supplemental information of the transform- 
ing microstructure, can therefore yield quantitative 
kinetic information. Computer modeling can also be 
used to devise and establish levels of confidence for 
methods that allow quick estimates to be made of 
important parameters for processing, such as the extent 
of the nucleation zone [26,27]. 

6. Conclusions and future directions 

A brief discussion of the polymorphic crystallization 
of glasses was presented. A new approximate expres- 
sion for finite-sample size effects on the transformation 
kinetics was discussed and was incorporated into a 
computer-based model for isothermal and nonisother- 


mal crystallization data. Good agreement between 
model predictions and experimental data the crystalliza- 
tion of lithium disilicate glass was demonstrated. 
Though a silicate glass was used for discussion, the 
general conclusions follow for all cases of polymorphic 
devitrification. Primary conclusions include the follow- 
ing: 

• the JMAK analysis is correct when applied within its 
range of validity; 

• for widely separated nucleation and growth regimes, 
common for silicate glasses, the Kissinger analysis of 
polymorphic crystallization data can provide an esti- 
mate of the activation energy for growth, but is 
invalid otherwise; 

• the Ozawa analysis gives reasonable estimates of the 
Avrami exponent for polymorphic transformations; 

• time-dependent nucleation can result in complicated 
devitrification behavior that cannot be understood 
with any simple method of analysis — realistic com- 
puter modeling is required; 

• a detailed study of the peak shape in DSC/DTA data 
can give more quantitative information than can be 
obtained from Kissinger-type methods; 
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m supplementary information of the transformation 

microstructure must be incorporated into any analy- 
sis of calorimetric data. 

The results presented here were for polymorphic 
devitrification by homogeneous nucleation. Clearly, 
they should be extended to the more common case of 
primary crystallization. Though diffusion-limited 
growth has been well studied and can be incorporated 
into the computer model in a rather straightforward 
manner, the correct treatment for the time-dependent 
nucleation of a cluster with a composition different 
from that of the parent phase is unclear. Here the 
embryo composition can change with cluster size and 
the chemical make-up of the regions surrounding the 
embryos will likely vary in a complicated manner, given 
the stochastic nature of cluster evolution. Often the 
classical theory can still be applied, changing only the 
Gibbs free energy [28,29] or the interfacial energy [30] 
with composition. When competition between diffusion 
in the parent phase and interfacial attachment arises, 
however, the classical theory is inadequate because the 
diffusive fluxes in the parent phase are coupled into the 
stochastic fluxes from the parent phase to clusters of 
the new phase. We have recently extended our model of 
time-dependent nucleation to treat this latter case. In- 
terestingly, the scaling of the nucleation rate and the 
transient time with the mobility as obtained from the 
classical theory changes. A more complex dependence 
on thermal history is therefore predicted, likely making 
the use of simple analytical expressions, such as the 
Kissinger analysis, even more invalid. 

Crystallization studies are rich, particularly when 
coupled with multi-step annealing treatments. More 
realistic models, such as the computer model described 
here, however, are required for the proper analysis of 
these data. The ready availability of the necessary 
computer processing power now makes this feasible. 
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